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0- " Abstract 

o: 

psj , The objective of the present paper is to investigate the solution of fully coupled 

mean-field forward-backward stochastic differential equations (FBSDEs in short) and 
to study the stochastic control problems of mean-field type as well as the mean-field 
stochastic game problems both in which state processes are described as FBSDEs. By 
combining classical FBSDEs methods introduced by Hu and Peng [Y. Hu, S. Peng, 



Solution of forward-backward stochastic differential equations, Probab. Theory Relat. 
Fields 103 (1995)] with specific arguments for fully coupled mean-field FBSDEs, we 
prove the existence and uniqueness of the solution to this kind of fully coupled mean- 
field FBSDEs under a certain "monotonicity" condition. Next, we are interested in 
optimal control problems for (fully coupled respectively) FBSDEs of mean-field type 
with a convex control domain. Note that the control problems are time inconsistent in 
the sense that the Bellman optimality principle does not hold. The stochastic maxi- 
mum principle (SMP) in integral form for mean-field controls, which is different from the 
classical one, is derived, specifying the necessary conditions for optimality. Sufficient 
conditions for the optimality of a control is also obtained under additional assump- 
tions. Then we are concerned the maximum principle for a new class of non-zero sum 
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stochastic differential games. This game system differs from the existing literature in 
the sense that the game systems here are characterized by (fully coupled respectively) 
FBSDEs in the mean-field framework. Our paper deduces necessary conditions as well 
as sufficient conditions in the form of maximum principle for open equilibrium point of 
this class of games respectively. 

Key words: Mean- field; forward-backward stochastic differential equation (FBSDEs); 
forward-backward stochastic control systems; stochastic maximum principle; non-zero sum 
stochastic differential game. 



1 Introduction 

In this paper, we consider the fully coupled forward-backward stochastic differential equa- 
tions (FBSDEs) of mean-field type 

X t = x+ [ E'[b(s,X' s X,K,X s ,Y s ,Z s )}ds+ [ E'[a(s,X' s X,Z' s ,X s ,Y s ,Z s )]dW s , 
Jo Jo 

Y t = $(X T ) + ^" E'[f(s,X' s X,Z' s ,Xs,Ys,Zs)}ds- J Z s dW s , (1) 



where b, f : Q x Q x [0, T] x R x R x R d x R x R x R d — ► R and a : Q x x [0, T] x 
RxRxR d xRxRxR d — ► R d satisfy the "monotone" condition introduced firstly 
by Hu and Peng [1] and W is a d— dimensional Brownian motion. Here, the coefficients 
E'[0(s, X' s , Yg, Z' s , X s , Y s , Z s )] (<fi = b, a, f), which are different from the classical coefficients 
of fully coupled FBSDEs, can be interpreted as 

E'[<p(s,X' s X,Z' s ,X s ,Y s ,Z s )] = [ 0(a/,w,a,X>O,^(wO,^(^) J ^-(w),i;(w),Z i (w))P(da/) 

Jn 

(Fully coupled) FBSDEs are encountered in the probabilistic interpretation (Feynman- 
Kac formula) of a large kind of second order quasi-linear PDEs, mathematical economics, 
mathematical finance and especially in the stochastic control problems (cf. [2]- [5]). There 
have been many results on the solvability of fully-coupled FBSDEs. Antonelli |6J first studied 
these equations, and he proved the existence and uniqueness with the help of the fixed point 
theorem when the time duration T is sufficiently small. Among others, to our knowledge, 
there exist three main methods to investigate the solvability of an FBSDEs on an arbitrarily 
prescribed time duration. The first one concerns a kind of "four step scheme" by Ma et al. 
[7| which can be regarded as a sort of combination of methods of PDEs and probability. The 
second one is the purely probabilistic method by Hu and Peng [lj, Peng and Wu [1], Yong [8j 
and Pardoux and Tang [9j. They required the "monotonicity" condition on the coefficients. 
The third one is motivated by the study of numerical methods for some linear FBSDEs (see 
Delarue [10J and Zhang [I lj). Delarue [10J relied on PDEs arguments, so its coefficients have 
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to be deterministic while Zhang [TT] imposed some assumptions on the derivatives of the 
coefficients instead of the monotonicity condition. 

Buckdahn, Djehiche, Li, and Peng [12] and Buckdahn and Li et al. [13] investigated a new 
kind of BSDEs-Mean-field BSDEs, inspired by Lasry and Lions [TJ]. In the present work, 
we adapt the methods developed by Hu and Peng [T] in order to establish the existence and 
uniqueness result for the fully coupled mean- field FBSDEs under the "monotone" condition. 
The two technical lemmas, aiming to prove the existence result of fully coupled mean-field 
FBSDEs, differ from the classical lemma in [1] because of the mean-field type. When the 
coefficients b, a and / do not depend on uj', the fully coupled equation ([I]) reduces to the 
standard one. So our result is nontrivially more general of PQ. 

We also consider stochastic optimal control problems and stochastic differential games 
(SDGs) in which the state variables are described by a system of mean-field FBSDEs. 
Mean-field control problems were recently studied by many researchers, such as Anders- 
son, Djehiche|15|. Buckdahn, Djehiche and Li [16], Meyer-Brandis, 0sendal, and Zhou [T7] 
and Li |18] . Andersson, Djehiche[15j use the methods in Bensoussan |19] to obtain the nec- 
essary conditions of the optimality of a control, i.e. they suppose that the control state 
space is convex so as to make a convex perturbation of the optimal control and obtain a 
maximum principle of local condition. Buckdahn, Djehiche and Li [16] get a Peng's type 
maximum principle for a general action space where the action space is not convex, using a 
spike variation of the optimal control. In Meyer-Brandis, 0sendal and Zhou [17] , a stochastic 
maximum principle of mean-field type in a similar setting is studied, but by using Malliavin 
calculus. Li [18J, also using the convex perturbation technology with the convex assumption 
for control domain, has a different controlled system and state equation of mean-field type 
from |15| . 

However, the results above are all on the forward control system. As far as we know, Peng 
|2U] originally studied one kind of forward-backward stochastic control system which has the 
economic background and could be used to study the recursive optimal control problem 
in the mathematical finance. He obtained the maximum principle for this kind of control 
system with the control domain being convex. Later, Shi and Wu [21] applied the spike 
variational technique to derive the maximum principle for fully coupled forward-backward 
stochastic control system in the global form and indicated that the control domain is not 
necessarily convex but the control variable can't enter into the diffusion term. In order 
to study the forward-backward stochastic control problem under the mean-field framework, 
we apply the convex perturbation methods introduced in Bensoussan [19J and analytical 
technique provided by jTS] to establish a necessary condition for optimality of the control in 
the form of the maximum principle for the (fully coupled respectively) mean-field forward- 
backward stochastic control system in which the state equation is mean-field FBSDE (fully 
coupled mean-field FBSDE respectively). The adjoint equation, playing an important role 
in deriving the SMP, is a (fully coupled respectively) mean-field backward SDE and has a 
unique adapted solution under the given assumptions with the help of the conclusion in [13] 
(or the conclusion in Theorem 3.1 respectively). Also, we obtain the corresponding sufficient 
condition, which can check whether the candidate optimal control is optimal or not. Our 
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results can be reviewed as an extension of Peng [20J and Li [18]. 

Inspired by Wang and Yu [22], which gave the maximum principle for non-zero sum 
differential games of BSDE system, we study the non-zero sum stochastic differential games 
(SDGs in short) of mean- field type. Differential games, originally studied by Isaacs [23], are 
ones in which the position, being controlled by players, evolves continuously. Fleming and 
Souganidis [24] were the first to study in a rigorous manner two-player zero sum SDGs. Their 
work has translated former results on differential games by Isaacs |23j . Friedman |25| . and, 
in particular, Evans and Souganidis |26j from the purely deterministic into the stochastic 
framework and has given an important impulse for the research in the theory of stochastic 
differential games. Next, the advances in SDGs appear over a large number of fields (cf. 

We notice that the game literature is mainly restricted to forward (stochastic) systems, 
i.e., these game systems are described by forward (stochastic) differential equations. Re- 
cently, Wang and Yu [22] concerned the theory of backward stochastic differential games 
and obtained the maximum principle as well as the verification theorem for non-zero sum 
SDGs of BSDEs in which game systems are described by BSDEs. It is remarkable that 
this topic about the forward-backward system is quite lacking in literature. To fill the gap, 
we investigate the theory of forward-backward SDG problems under the mean-field frame- 
work. Similar to our stochastic control problems, we study the SDGs with the state equation 
having two different forms: mean-field FBSDEs and fully coupled mean-field FBSDE. By 
virtue of an argument of the convex perturbation, we deduce the stochastic maximum prin- 
ciple for the equilibrium point of Problem (FBNZ) (Problem (CFBNZ) respectively), which 
gives the candidate equilibrium points. By extending classical approaches to the mean-field 
framework, we prove, under some restrictive assumptions (but comparable with those in 
the classical case), the sufficiency of the necessary conditions. It is necessary to point that 
our SDGs conclusion not only extends the result of Wang and Yu [22] but also includes 
the situation where the state equation of the stochastic game system is classical (i.e. in no 
mean- field form) FBSDE (fully coupled FBSDE respectively). 

Our paper is organized as follows. Section 2 recalls some elements of the theory of 
FBSDEs and mean-field BSDEs which are needed in what follows. Section 3 investigates 
the uniqueness and existence of the solution of fully coupled mean-field FBSDEs under the 
"monotonicity" condition in which two technical lemmas are used to prove the existence 
result. In Section 4, we study the forward-backward stochastic control system of mean-field 
type. Specifically, the maximum principle, specifying the necessary condition for optimality, 
is deduced and we get, under additional assumptions, the corresponding sufficient condition 
which can check whether the candidate control is optimal or not. Similar results about 
fully coupled forward-backward stochastic control system of mean-field type are obtained in 
Section 5. Following the idea introduced in Section 4 and Section 5, we analyze the non- 
zero sum stochastic differential games of FBSDEs and fully coupled FBSDEs in Section 6 
and Section 7, respectively, and derive the necessary condition in the form of the maximum 
principle as well as the sufficient condition- verification theorem for the equilibrium point. 
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2 Preliminaries 



Let (Q, F, Ft, P) be a given complete filtered probability space on which a (/-dimensional 
standard Brownian motion W = (Wt)t>o is defined. By F = {Ft,0 < t < T} we denote the 
natural filtration of W augmented by P— null sets of F, i.e., 

F t = a{W s ,s<t}VAf P , te[0,T), 

where Hp is the set of all P-null sets and T > is a fixed time horizon. 

We shall introduce the following two processes which can be used frequently in what 
follows: 



«Sf(0, T; R) : = 

(0 4 )o<t<T real-valued F — adapted cadlag process : E sup \(j) t 

\-a<t<T 



< +oo 







" f\<Pt?dt 


< +oo| 







Mj(0, T; R n ) : = j OJ <t<T K n - valued F - adapted process : E 

2.1 The classical FBSDEs 

We first recall some results on FBSDEs, for its proof the reader is referred to Hu and Peng 
PQ. The FBSDEs they considered has the form 

x t = x + f*b(s,x s ,y s ,z s )ds + f* a(s,x s ,y s , z s )dW s , 

Vt = g{x T ) + J t T f(s, x s , y s , z s )ds - z s dW s , t G [0, T\. 

Function b, f : x [0, T] x R x R x R d ->■ R, a : Vt x [0, T] x R x R x M d ->• M d with 
the property that b(t, x, y, z) te[0:T] ,a(t, x, y, 2)te[o,T] and f(t, x, y, z) te [o,T] are F-progressively 
measurable for each (x, y, z) G R x R x M d . 

Some notations and conditions are needed before giving the existence and uniqueness 
of the solution of such FBSDEs. Let <, > denote the usual inner product in R n , and for 
u = (x, y, z) G R x R x R d , we define 

F(t,u) := (-f{t,u),b(t,u),a(t,u)). 

(HI) (i) For each u = (x, y, z) G R x R x R d , F(-,u) G .M 2 (0, T;lxlxR d ), and for each 
x G R, g(x) G L 2 (f2, J-y; K); there exists a constant C\ > 0, such that 

\F(t, ui) — F(t, w 2 )| < Ci|«i — M2I, P — a.s., a.e. £ G M + , 
VUi G I x 1 x K d (i = 1,2) ; 
|g(xi) — 5 l (x2) I < ci|xi — x 2 \, P — a.s., V(xi, x 2 ) G R x R. 
(ii) There exists a constant c 2 > 0, such that 

< F(t,ui) - F(t,u 2 ),ui - u 2 > < —c 2 \ui - u 2 \ 2 , 

P - a.s., a.e. t G R + , V u { G E x R x ]R d (« = 1, 2) ; 

< — g(x 2 ), x\ — x 2 > > c 2 \xi — x 2 \ 2 , P — a.s., V (x±, x 2 ) Glxl. 
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Lemma 1. Let assumptions (HI) hold, then there exists a unique adapted solution (x,y,z) 
for the FBSDEs (1 ) 

2.2 Mean-field BSDEs and McKean-Vlasov SDEs 

This section is devoted to the recall of some basic results on a new type of BSDEs, the 
so called mean-field BSDEs; the reader interested in more details is referred to Buckdahn, 
Djehiche, Li, and Peng [12] and Buckdahn and Li et al. [13J. 

Let (n,J,P) = (Q x f2, P (g) P, P (g) P) be the (non-completed) product of (fl, P, P) with 
itself. We endow this product space with the filtration F = {P t = P £g> P(,0 < t < T}. 
Any random variable £ G L°($7, J 7 , P) originally defined on Q is extended canonically to Q : 
g{uj',uj) = G fi = O x fi. For any G L 1 ^, P,P) the variable 0(-,w) : ft R 

belongs to L 1 (fi, J 7 , P) , P(du) — a.s.; we denote its expectation by 

E'[0(.,w)] = / 9(u',uj)P(duj'). 
Jn 

Notice that E'[0] = E' 1 ^, J",P), and 

P[6>] = [ 6dP= [ E'[e(-,co)}P(duj) = E[E'[0]]. 
Jn in 

The driver of mean-field BSDE is a function / = /(a;', u, t, y, z,y, z) : O X [0, T] x 1 x l* 1 x 
which is F-progressively measurable for all (y, z, y, z), and satisfies the following 

assumptions. 

(H2) (i) There exists a constant C > such that, P-a.s., for all t G [0, T], 2/1, 2/2, 2/1, 2/2 G R, 
21,^,^1,^2 e |/(*, 2/1, zi, 2/1, zi) -f(t, y 2 , Z2,V2,z 2 )\ < C(|j/i -2/2I + 1^1-^1 + 
2/2 1 + \zi - z 2 \). 

(ii) /(.,0,0,0,0)e^ f (0,T;R). 
The main result about mean-field BSDEs of Buckdahn and Li et al. [13] is: 

Lemma 2. Under the assumptions (H2), for any random variable £ G L 2 (Q , Tt , P) , £/ie 
mean- field BSDEs 

F * = £ + / W\f(sX,Z'„Ys,Zs)]ds- J Z s dW s , 0<t<T, (2) 

/ias a unique adapted solution 

(Y t ,Z t ) e S 2 (0,T;R) x M 2 F (0,T;R d ). 
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Remark 3. The driving coefficient of (0|) has to be interpreted as follows: 

E'[f(sX,Z' s ,Y s ,Z s )]^) 
= E'[f{s, Y s (u'), Z s (u'),Y s (lu), Z s (u))} 

f(u', u, s, Y s (u)'), Z s (u'), Y s (u), Z s (u))P(du'). 



We shall also consider McKean-Vlasov SDEs (see, e.g., Buckdahn and Li et al. [13]) ■ Let 



b : tt x [0, T] x R x R d -> R and a : ft x [0, T] x 
supposed to satisfy the following conditions: 



be two measurable functions 



(H3) (i) b(-,x,x) and a(-,x,x) are F-progressively measurable continuous processes for all 
1,16 8, and there exists some constant C > such that 

\b(t,x,x)\ + \a(t,x,x)\ < C(l + \x\ + \x\), a.s., 

for all < t < T; 

(ii) b and a are Lipschitz in x, x, i.e., there is some constant C > such that 

\b(t, Xi,xi) - b(t,x 2 ,x 2 )\ + \(r(t, xi, xi) - a(t,x 2 ,x 2 )\ < C(\xi - x 2 \ + \x\ - x 2 \),a.s. 

for all < t < T, X\, x 2 ,xi, x 2 G M. 

The McKean-Vlasov SDEs parameterized by the initial condition (t, G [0, T]xL 2 (Q, J- t , P; M.) 
is given as follows: 



dXl't = E' 



1 s 



ds + E' 



We recall that, due to our notational convention, 

E' [ft (a, (X^)' , X*'^] (w) = J b (to', co, s, X 1 / {to') , X 1 / (w)) P(dt/), w G fi. 

Lemma 4. Under Assumption (H3), SDEs () has a unique strong solution. 

Remark 5. From standard arguments we also get that, for any p > 2, there exists C p G R, 
which only depends on the Lipschitz and the growth constants of b and a, such that for all 
t G [0,T] and(,(' p (n,F t ,P;M), 



E 
E 
E 



sup 

L t<s<T 

sup 

t<s<T 



\ x t,C _ x t /\P\F t 



<c p \c-C p l 



a.s. 



\X\' Q \ P \Ft 



<C p (l + |Cn, as., 



(3) 



sup ix^-ci^] <c p (i + \c\nsi 

L t<s<t+S J 



P-a.s., for all 5 > with t + 5 <T . 

These, in the classical case, well-known standard estimates can be consulted, for instance, 
in Ikeda and Watanabe [30] (pp. 166-168) and also in Karatzas and Shreve f31f(pp. 289-290). 
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3 Fully coupled Mean-field FBSDEs 



In this section, we shall investigate a new type of FBSDEs called fully coupled mean-field 
FBSDEs as follows: 

X t = X + £ E'[b(s, X' s , Y>, Z' s , X s , Y s , Z s )]ds + f* E'[a(s, X' s , Y>, Z' s , X s , Y s , Z s )]dW s , 

Y t = $(x T ) + f t T E'[f( s , x'„ y;, z' s , x s , y , z s )]ds - f t T z s dw s , t e [o, r]. 

Here the processes X, Y, Z take values in R, R, R d respectively; and b, a, $ and / take values 
in R, R d , R and R respectively. 

Remark 6. The driving coefficient here has to the same interpretation as Lemma 2: 

E'[^{ S ,X' s X,Z' s ,X s ,Y s ,Z s )]{oj) 
= E'[^(s, X s (u/),Y 8 (u/), Z s (u/),X 8 (u>), Y s (oo), Z s {u))] 

= ( iP(uj', oo, s, X s (uj'),Y s (uj'),Z s (uj'),X s (uj), Y s {oj),Z s {u))P{dJ). 
Jn 

for ip = b,a, f. 

For convenience, we will use the following notations in this section: Let <, > denote 
the usual inner product in R n and we use the usual Euclidean norm in R n . For = 

(x, y, z, x, y, z) G R x R x R d x R x R x R d , 

F(t,e) = (-f(t,e),b(t,e),a(t,&)). 

Now we give the standard assumptions on the coefficients of mean-field FBSDE: 

(H4) For each 6GRxRxR d xRxRxR d , F(-, 0) G M 2 (0, T;RxRxR d xRxRx R d ), 
and for each x G R, g(x) G L 2 (Vl, J 7 , R); there exists a constant C > 0, such that: 

\F(t, 0i) - F(t, 2 )| < C|0i - 2 ], P - a.s., a.e. t G R+, 

©i = (x h y t , z h xt, yt, Zi) G R x R x R d x R x R x R d , (i = 1,2), 

and 

|$(xi) -$(x 2 )| < C|xi -ar 2 |, P - a.s.,V(x 1 ,x 2 ) G R x R. 

The following monotone conditions are our main assumptions: 

(H5) For Qi = (x h y h z h x^y^Zi) elxRxt rf xlxEx R d , let m = (x i: y t , z { ) G R x R x R d , 
then 0j = (ui,Ui) (i = 1,2). We assume that 

E < F(t, 6i) - F(t, e 2 ),tii - m 2 > < -CiEflui - m 2 | 2 ) P - a.s., a.e. i G R + , 
< $(x 1 )-$(x 2 ),xi-X2> > ^ 1 \xi-x 2 \ 2 , P - a.s., V(x 1 ,x 2 ) G R x R, 

where Ci and /i 1 are given positive constants. 



8 



For the mean-field FBSDE (j3J), we have the the following main result of this section. 

Theorem 7. Under the assumptions (H4) and (H5), there exists a unique adapted solution 
(X, Y,Z) for mean-field FBSDEs (TJj. 

The proof of this theorem is similar to that of Theorem 3.1 in [I] except the mean-field 
term. However, to be self-contained, we intend to give the proof. Before giving the proof 
of this theorem, we need the two technical lemmas below whose proof will be given in the 
sequel. 

Lemma 8. Suppose that (^(-), </>(■), ip(-)) G ^ 2 (0,T;lx R d xl), £ G L 2 (Q,J" T ;M), then 
the following linear mean- field forward- backward stochastic differential equations 



(4) 



X t = X Q + f ( - E'[Y S '] -Y. + 7 (a))ds + f ( - E'[Z' S ] -Z s + <j>{s))dW s , 
Jo Jo 

Y t = C + X T + [ \E'[X' a ]+X a -<p(s)]d8- I Z s dW s , (5) 
Jt Jt 

have a unique adapted solution: (X, Y, Z) eM 2 (0,T;lxlx R d ). 
Now, we define, for any given aGl, 

b a (t, x, y, z, x, y, z) = ab(t, x, y, z, x, y, z) + (1 - a)(-y - y), 
a a {t, x, y, z, x, y, z) = aa(t, x, y, z, x, y, z) + (1 - a)(-z - z), 
f a (t, x, y, z, x, y, z) = af(t, x, y, z, x, y, z) + (a - l)(-£ - x), 
=a$(x) + (l-ai)(x). 

and consider the following equations: 

X t = X + f [b a (s,A s ) + j(s)]ds+ [ [a a (s,A s ) + ( f ) (s)}dW s , 
Jo Jo 

Y t = ($ Q (X T ) + e)+ / [f a (s,A s )-<p{s)]d S - I Z s dW s , (7) 

Jt Jt 

where we use the notation 

A* = (X s , Yg, Z s , X s , Yg, Z 8 ), 

and 

$(s, A s ) = E'[i;{X'g, Yg 1 , Z' s , Xg, Yg, Z s )l 
for ip = b,o~, f. Then we can rewrite 

b a (s,A s ) = ab(s,Ag) + {l-a)(-E'[Y s '}-Yg), 
a a (s,A s ) = aa(s,Ag) + (l-a)(-E'[Z' s }-Zg), 
r( Sl A s ) = af(s,A s ) + ( a -l)(-E'[X' s ]-X s ). 



(6) 
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Lemma 9. For a given a G [0, 1) and for any (■),(/)(■), <£>(•)) G M 2 (0, TjlxR^xl), £ e 
L 2 (Q, Tt, P; R) ; assume that Eqs (0|) and ([?)) nave an adapted solution. Then there exists a 
5q G (0,1) which depends only on Cj, C2 and T ; snca iaai /or a// a G [«oi a o + ^o]> an d f or 
any (/>(■), <p(-)) G A1 2 (0,T;R x R d xR), £ G L 2 (fi, J" T , P; R), Pgs ^ and (TJ] nave an 

adapted solution. 

Proof of Theorem 7. 

Uniqueness. If U = (X,Y,Z) and {7 = (X,Y,Z) are two adapted solutions of ((3]), we 

set 

(X\ Y', Z', X, Y, Z) = (X' - X', Y' - Y' , Z' - Z\ X - X, Y - Y, Z - Z), 
b(t) = b(t,U',U) - b(t,U',U), 
a(t) = a(t,U',U) - a(t,U',U), 
f(t) = f(t,U',U)-f(t,U',U). 

From Assumption (Al), it follows that {X t } and {Y t } are continuous, and 

E( sup \X t \ 2 + sup \Y t \ 2 ) < +00. 
te[o,T] te[o,T] 

Applying the Ito's formula to X t Y t on [0, T], we have 

E[($(X T )-Q(X T ))X T ] 

= E^ {E'lb(t)]Y t -E'lf(t)}X t + E'la(t)]Z t }dt. 

= eJ^ { < {-E'[f(t))Mkt))MHt)]),{XtXZt) > }dt 
By assumptions (Al) and (A2), we get then 

fi 2 \X T - X T \ 2 < E[($(X T ) - $(X T ))X T ] < -dE [ \U- U\ 2 dt. 

Jo 

So, we get U = U. 

Existence. According to Lemma 8, we see immediately that, when a = 0, for any 

(7(0^(0, 

¥>(■)) G M 2 (0, T;lxl d xl),(G L 2 (Q, T T , P; R), Eqs © and © have an adapted solution. 
From Lemma 3.2, for any (7O), <£(•), tp(-)) G A1 2 (0,T;R x R ' x I), ( G L 2 (fi, J" T , P; R), we 
can solve Eqs (jBJ) and (j7|) successively for the case a G [0, o" ], [o" , 25 ], • • • ■ When a — 1, for 
any ( 7 (-), <£(•), </?(•)) G A4 2 (0,T;R x R d x R), f e L 2 (£7, J^, P; R), the adapted solution of 
Eqs. ([6]) and (J7J) exists, then we deduce immediately that the adapted solution of Eqs. ([3|) 
exists. □ 
Proof of Lemma 8 
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Proof. We consider the following BSDEs: 

K = Z + J [-l'R]-y,-^)+7(#- J (E'[Z' a ] + 2Z a -<f>(s))dW a . 

By Lemma 1, the above equation has a unique adapted solution (Y, Z). 
Then we solve the following forward equation 

X t = x + f ( - E'[X' S ] -X s - E'K'] - Y s + 7 (s))da + / ( - E'[Z' S \ -Z s + <j>{s))dW s 
Jo Jo 

and set Y = Y + X, Z = Z, we get 

X* = x+ f (-E'[Y^]-Y s + 1 (s))ds+ [ (-E'[Z' S ]- Z s + <f){s))dW s , 
Jo Jo 

t + j [E'[X' S ] - E'K'] + X s - Y s - ^(a) + 7 ( fl )]d fl 

{E'[Z' s \+2Z B -(j){s))dW s , 
J ( - E'K] -Y s + j{s))ds + J (-E'[Z' S ]- Z s + <j)(s))dW s . 



Y t -X t 



Xt — Xt 
Then we have 

Y t = i + X T + J ^[X' s ]+X s -v{s)\ds- j Z s dW s . 

So (X, Y, Z) is a solution of Eqs. (jl]) and flH]). Thus the existence is proved. 

As for uniqueness, it only has to use the method of the proof of uniqueness in Theorem 
3.1 and we omit it. □ 

Proof of Lemma 9 

Proof. For simplicity, we set 

XT = {X\Y\Z l ), 

A° = ((X y,{Y )',{Z y,X ,Y ,Z ) =0, 

A i = ((x i y,(Y i y,(z i y,x i ,Y i ,z i ), 
x i+1 = - x i = ((x i+1 y, (Y t+1 y, (z i+1 y, x i+ \ y i+1 , z i+l ), 

= ^ _ Jji = y i+ l^ + 

for all ieN+. 
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For any given a G [0, 1] and any 5 > 0, we solve iteratively the following equations: 
+ jf (a<*°(s, Al +1 ) +5[Zl + E'[(Zl)'] + a(s, A*)] + 0( S ))dW s , 



(8) 



+S[f(s,A i s )-n(Xiy}-Xl}- ( p(s))ds- J" Zl +l dW s . (9) 



Applying the Ito formula to X l t + Y^ , on and noticing that E'[Y'] = E[Y], we have 

E(($ ao (A^ +1 ) - $ ao (X^))A^ +1 ) 
= -<JE[($(X£) - §{X l T l ) - X l T )X l T +1 ] 

+E £ 1 1>; +1 [6«o ( a , A; +1 ) - b a » (s, A*)] - If 1 [r° ( S , Al +1 )- r(s, Al)] 
+Z'+ l [a«°(s, Al +l ) - a«°( S , A^ds + SE £ [*? + E[*?] + &(«, A*j - 6( a , A^ 1 ) 

[/(*, A J j - /(,, Ar x ) - X* - EK]] + Zi +t [Zl + E[Zi\ + a(s, Al) - a(s, A*" 1 )] \ds 

= -6E[($(X l T ) -^(X^ 1 ) - X l T )X l T +1 ) +eJ^ (^F a °(s,X l s +l ) - F a °(s, Al),Ul +l ^ ds 
+5E £ (Ul + E[Ul] + F(s, A*) - F(s, A^ 1 ), C/f 1 } ds. 
Using the definition of F(t,A), we have the following inequality: 



F( fl ,Ai)-F(a,AjT l ) = E'[F( S ,A;)-F( S ,Ar^ 

< CE'IA^-A^ 1 ! 

= C(E[|f> i |] + |f> i |). (10) 
From assumptions (Al), (A2), Eq (fTUj) and the notation 

F a °(s, A* +1 ) = a F(s, Al +1 ) - (1 - a ) (E[(f>; +1 )'] + 
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for any a £ [0, 1], we deduce easily that 

(frao + 1 - a )E[X l T +1 } 2 + (dao + 2- 2a )E [ \Ut +1 \ 2 dt 

Jo 

< 5E\X % T X l + l ] - 5E[($(X l T ) - ^(Xif 1 ))^ 1 ] +5E f < U l t , U' t +l > dt 

Jo 

+5E [ < E[U l t ],Ui +1 >dt + 5E [ < F(t, Aj) - F(t, A*" 1 ), U l t +1 > dt 
Jo Jo 

< 25(1 + C) (e\X 1 t \\X 1 t +1 \ + E^ \U l t \\U l t +1 \d?j . 
In virtue of 

min (/i 1 «o + 1 — &o, Ci^o + 2 — 2a ) > A* 2 = min (1, Ci) > 0, 



the above inequality yields 



E[X 1 t +1 ] 2 + e[ \U l t +1 \ 2 dt < 2 ^ 1 + C W |X^||Xj+ 1 |+ / \Ui\\Ui +1 \dt 

Jo ^2 V Jo 

For e = 4<5( ^ c) > 0, with the help of ab < + £ 6 2 , we can derive 

E[l^ +1 ] 2 + E f \Ut +1 \ 2 dt< 
Jo 



25(1 + C) i E ( |^ )2+ / |^|2^j 



^2 

Note that there exists a constant C > which depends only on C and T, such that 

E|A^| 2 < C'eJ^ (\U l s \ 2 + lUi'^ds, Vi > 1. (12) 

Indeed, for i > I, 

e\x* t \ 2 = mj |x*|{^( s ,a*)-^>(^ 

{^( S ,a:) -o-o^Aj- 1 ) +5(ir 1 +E / [(ir 1 )i +^,Ar 1 ) -^oar 2 ))} 

< E^ \X l s \ 2 dt + 2E^ {(b a "(s, Al) -^(s,^- 1 )) 2 + {a ao (s,XD - a a « (s , A^ 1 )) 2 } ds 

{y;- 1 + E'^y;- 1 )'] + l( s , a*- 1 ) - 5(a, a;- 2 ) } ds 

+25 2 E £ {zr 1 + E'[(Zlr 1 )'} + A^" 1 ) - a(s, A\7 2 )} 2 ds 

< CE f T {\Ul\ 2 + |L>r 1 | 2 + E[|l > ;- 1 | 2 + \Z\- l \ 2 ] +E(|^| 2 )+E(|L>r 1 | 2 )}ds. 



2 

ds 
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Similar to (|10|) . in the last inequality, we use the fact that 

feArV&^Ar 2 )! 2 < E' ( [6(s, A* -1 ) - 6(s, A* -2 )] 2 ) 

< c(E(|t>r i i 2 ) + i^- i i 2 ), 

as well as 

^(s^D-^^Ar 1 ) < a C|A i s -Ar 1 | + (l-a )(-E[y l ]-n ) 

where c is a constant which depends on C. One can show that a(s, Ag _1 ) — a(s,A* -2 ) and 
a a °(s,A*) — <7 ao (s, A* -1 ) have the similar results. By a standard method of estimation, the 
desired result (I12p can be derived easily. 

From (ITT]) and ( I12p , we know that there exists a constant K > which depends only on 
C, Ci, /x 1 and T, such that 

E^ |t> s m | 2 (is < K5 2 (e {|L>f + |^rT}^ 

Hence there exists a <5q G (0, 1), which depends only on C, Ci, ^ and T, such that when 
< 5 < 6 , 

E f \Ui +1 \ 2 ds < ^E [ \Ut\ 2 ds + lE [ lUjr^ds. 



ei [\ut +1 \ 2 + ^-\ui\ 2 )ds< -e I (if>:i 2 + -|i>r i i 2 ) ( i,s. 



That is 



Repeat the above inequality as many times as you desire, there holds 
"T , i , /i \ »-i rT 



E 



I (\u: +1 \ 2 + \m 2 )ds<(iy eJ q (i^i 2 +i|^i a )ds, i>i. 



It turns out that [/' is a Cauchy sequence in A4 2 (0, T;Mxlx M. d ) and its limit is denoted 
by U = (X, Y, Z). Passing to the limit in Eqs.® and we see that, when < d < Sq, U = 
(X, Y, Z) solves Eqs. ([6]) and ([7]) for a = (Xq + S. The proof is completed. □ 

The condition (A2) can be replaced by the following condition. 

(H6) For Qi = {x h y h z h x h y h Zi) E Rxtxl^xlxlx M. d , let m = {x h y h Zi) E lxRxR d , 
then 0j = (v,i,Ui) (i = 1,2). We assume that 

E < F(t, ©i) — F(t, Q 2 ), u± — u 2 > > CiE(|«i -u 2 \ 2 ), P — a.s., a.e. t E IR + , 

< <b(xi) — §{x 2 ),x\ — x 2 > < - fi^xi - x 2 \ 2 , P - a.s., V (xi, x 2 ) E R x R, 

where Ci and // x are given positive constants. 
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We have another parallel existence and uniqueness theorem for mean-field FBSDEs. 

Theorem 10. Let (H4) and (H6) hold. Then there exists a unique adapted solution (X,Y,Z) 
of mean- field FBSDEs (TJj. 

The method to prove the existence is similar to Theorem 7. We now consider the following 
(USD for each a G [0, 1] : 

dX* = [ab(s, A s ) + j(s)]ds + [aa(s, A s ) + cj)(s)]dW s , 
-dY s = [ - (1 - a)c 2 X s + af(s, A s ) + cp(s)] ds - Z s dW s , (13) 
X a = a, Y« = a$(X T ) - (1 - a)X T + £, 

where 7, (j) and ip are given processes in M 2 (0,T) with values in E, R d , and E, resp., 
£ G L 2 (Q, JFt, P)- Clearly the existence of (TTBl for a = 1 implies the existence of FBSDEs 
([3]). From the existence and uniqueness of SDEs and BSDEs, when a = 0, the equation (TTBl 
has a unique solution. 

In order to obtain this conclusion, we also need the following lemma. This lemma gives 
a priori estimate for the existence interval of (fT3|) with respect to a G [0, 1]. 

Lemma 11. We assume (H4) and (H6). Then there exists a positive constant 8q such that 
if, a prior, for a a® G [0, 1) there exists a triple of solution (X a °, Y a °, Z a °) of < f73j) . then for 
each 5 G [0, 5 ] there exists a solution (X ao+s , Y ao+s , Z aa+S ) of FBSDEs (Ef fora = a + S. 

Proof. We use the notations 



u = 




u = (x,y,z), 


U = 


(X,Y,Z), 


U=(X,Y,Z), 


= 


(x' : y',z',x,y,z), 


6 = (x\y',z,x,y,z), 


A = 


(X',Y',Z',X,Y,Z), 


A = (X>,Y',Z',X,Y,Z), 


e = 


e-e, 


A = A- A, 


u = 


(x,y,z) = {x-x,y- 


V,z- z), 


U = 


(X,Y,Z) = (X-X, 


Y-Y,Z-Z). 



Since for (7, <p, if) G M 2 {0,T;R x R d x E), £ G L 2 (^,7" T ,P), a G [0,1) there exists a 
unique solution of (lT3"j) . thus, for each Xt G L 2 (f2, J-^, P) and a triple it s = (x s ,y s ,z s ) G 
M 2 {0, T; E x M ' x E) there exists a unique triple £/, = (X s , F s , Z a ) G .M 2 (0, Tjlxf'xl) 
satisfying the following FBSDE: 

dX s = [a b(s, A s ) + 5b(s, Q s ) + j(s)] ds + [a a(s, A s ) + Sa(s, Q s ) + 0(s)] dW s , 
-dY s = [ - (1 - a )dX s + oof (a, A a ) + 5{C lX . + f(s, Q s )) + f(s)]ds - Z s dW s (U) 
X = a, Y T = a $(X T ) + (a - l)X T + 5(®{x T ) + x T ) + f . 
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We now proceed to prove that, if 5 is sufficiently small, the mapping defined by 

I ao +s(u x x T ) = U x X T : 
M 2 (0, T;lxR d xl)x L 2 (VL, JF t , P) ->■ M 2 (0, T;RxR d xR) x L 2 (VL, JF t , P) 

is a contraction. 

Let u = (x, y, z) E M 2 (0, T;RxR d x R), and let U x X T = I ao+s (u x x T ). 
Using Ito's formula to X S Y S yields 

a E < $(X T ) - <5>{X T ),X T > +{a - l)E\X T \ 2 + 5E < $(x T ) - $(x T ) + x T ,X T > 

= -E [ Xs\a [f(s,As)-f(s,As)] + S[f(s,e s )-f(s,e s )]-(l-a )C 1 X s + 6C 1 x s 
Jo 1 J 

/ Y s \a [b(s, A s ) - b(s, A s )] + 5[b(s, Q s ) - b(s, Q s )] 
Jo l 



ds 



+E 



ds 



+E 



ds 



= a E 



a [a(s, A s ) - a(s, A s )} + 5[a(s, <d s ) - a(s, Q s )] 

,[ < F(s,A s ) -F(s,A s ),U s > ds + (l -a )CiE / \X s \ 2 ds - 5dE [ < X s 
Jo Jo Jo 

+SE [ < F(s, Q s ) - F(s, Q S ),U S > ds. 
Jo 

From (H4) and (H6), we can get 

(frao + (f - ao))E[|X T | 2 ] + CiE / \X s \ 2 ds + da Q E [ (\Y S \ 2 + \Z s \ 2 )ds 

Jo Jo 

< 5E< $(x T ) - $(x T ) + x T , X T > +5CiE f < X s , x s > ds 

Jo 



x s > ds 



-SE / < F(s, Q s ) - F(s, G 8 ), U s > ds 
Jo 

< d^Edxrl 2 + \X T \ 2 ) + 5K 1 E (j 
This means 



Us + 



\Us\ 2 ) 



ds. 



fj,E[\X T \ 2 ] + CiE / \X s \ 2 ds < 5K{K(\ 
Jo 



x T \ 2 + \X T \ 2 ) + SK ± E 



u s \ 2 + 



\U s \ 2 )ds, 



where /^cko + (1 — «o) > A 4 = min(l,/i 1 ) > 0. 

On the other hand, for the difference of the solutions (Y, Z) — (Y — Y, Z — Z), we apply 
the usual technique to the BSDE part: 



E (\Y S \ 2 + \Z s \ 2 )ds < K 2 5E / \u s \ 2 ds + K 2 5E\x T \ 2 + K 2 E \X s \ 2 ds + C\E\X T 



f 

Jo 
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Here the constant K% depends on the Lipschitz constants C, Ci, and T. 
Combining the above two estimates, it is clear that, we have 

E[\X T \ 2 }+E J \U s \ 2 ds<5KE^J^ \ii s \ 2 ds+\x T \ 2 y 

Here the constant K depends only on C. C\, fj, 1 , Ki, and T. We now choose So = . It is 
clear that, for each fixed 5 G [0, 5q], the mapping I ao +s is a contraction in the sense that 

E[\X T \ 2 ]+E [ \U s \ 2 ds < Ue [ \u s \ 2 ds + E\x T \ 2 ). 

Jo 2 V J / 

It follows that this mapping has a unique fixed point U ao+s = (X ao+s , Y ao+s , Z ao+s ) which 
is the solution of ([13]) for a = a + 5. The proof is complete. □ 

We now give the proof of Theorem 10. 

Proof of Theorem 10. The uniqueness is obvious from Theorem 7. When a = 0, the 
equation (Tl3|) has a unique solution. It then follows from Lemma 3.3 that there exists a 
positive constant Sq depending on Lipschitz constants C, C\, K\ and T such that, for 
each 5 G [0,(5>o], equation (TT3"|) for a = + 5 has a unique solution. We can repeat this 
process for iV-times with 1 < N5q < 1 + 6q. It then follows that, in particular, FBSDEs (1131) 
for a = 1 with £ = has a unique solution. The proof is complete. □ 

Theorem 10 can ensure the existence and uniqueness of solution to the adjoint forward- 
backward systems in Section 5 and Section 7. 



4 Stochastic maximum principle in mean-field controls of 
FBSDEs 

In this section, we study the stochastic maximum principle for mean-field control problem 
of FBSDEs. The action space U is a non-empty, closed and convex subset of R k (k G N + ), 
and we define the admissible control set as 

IA = {v t G Zv|(0,T; U)\vt(u',u) : [0, T] x Q x Q U, t G [0,T]}. 

For any v(-) G U, we consider the following forward-backward stochastic control system of 
Mean- field type: 

dX t = W[b(t, X;, X u v t )\dt + E'[a(t, X>, X t , v t )}dW t , 
X(0) = Xq, 

-dY t = E'[f(t, X' t , y/, Z' t , X t , Y u Z u v t )]dt - Z t dW t , 
Y T = $(X T ), 
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where 

b : [0,T] x R x R x U -»• R, 

a : [0, T] x R x R x U -> R d , 

/ : [0, T] x R x R x R d x R x R x R d x U -»■ R, 

$ : R — > R. 

The optimal control problem is to minimize the following expected cost functional over U: 

fT 



J(«(0) = E(7 X t ', y/, Z' t , X t , Y t , Z t , v(t))]dt 

+E^(X r )+ T (y(0))), (15) 



where 



R — )■ R, 
R -> R, 

[0, T] x R x R x R fc x R x R x R fc x U -> 



An admissible control u ElA is said to be optimal if 

J(ii) = min J(v) . (16) 

Now we give the following conditions in this section. 

(Al) The given functions b(t, x, x, v), a(t, x, x, v), f(t, x, y, z, x, y, z, v), h(t, x, y, z, x, y, z, v),g(x) 
and j(y) are continuously different iable with respect to all of their components respec- 
tively. 

(A2) All the derivatives in (Al) are Lipschitz continuous and bounded. 

For any admissible controls v(-) G U, due to Lemma 2, the mean-field FBSDEs (j3j) admits 
a unique solution under assumptions (Al) and (A2), which is denoted by (Xt, Y t , Z t ). 



4.1 Variational equations and variational inequality 

Let u(-) be an optimal control and (X u (-), Y u (-), Z u (-)) be the corresponding state trajectory 
of stochastic control system. For any < 9 < 1, we denote by (X® ,Y t e , Z®) the state 
trajectory corresponding the following perturbation u 9 t of u t . 

u\ = u t + 9(v t -u t ), v t eU. 
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Since U is convex, then u (•) is also in U. Let (&(•), m(-), n(-)) be a solution of the variational 
equation 



dh = E' 



6 5 (t, (X t ")*> t )(^)' + W)?S«t)*t 



tft + E' 



dW t 



(17) 



+<r ie (t, (X«)f >Ut )fci + a v (t, (X?)',X t ,u t )(v t - ut) 
k = 0, 

dm t = -E'[f £ (t,u)(kt)' +_f x {t,u)k t + fy(t,u)(m t y + f y (t,u)m t + f s (t,u)(n t y 

+f z (t,u)n t + f v (t,u)(v t - ut)]dt + n t dW t , 
m T = k T $ x (X T ). 



where we use the notation f(t,u) = f(t, (X™] 
Set 



iu\iu\lu vu yu „, \ 
t It It i Y t > Z t > U tl 



Zf 



X?) - k 



t- 



e-\xl 

e-\Y t 6 -Y t u )-m u 
e-\z e t - Z?) - n t . 



Then, we have the following convergence result: 
Lemma 12. We suppose (A3) and (A4) hold. Then 

12 = 0, 

2 = o, 

| 2 = 0. 



lim sup E\Xt 

e ^°0<t<T 

lim sup E\Y t 9 

0->O(Kt<T 



lim sup E\Z t 

0<t<T 



(19) 



(20) 



Proof. Since the coefficients in linear mean-field FBSDE (|17|) and (118]) are bounded, it follows 
from Proposition 1.2 in [13J that there exists a unique solution (k(t), m(t), n(t)) for equations 
([IT]) and (USD. 

The proof for the convergence of X 6 t can be found in Lemma 3.2 of [18J. We need only 
to deal with Y t e and Zf. From the definition of Y t e , it fulfills the following BSDE, 



-dYf = --{dYf - dY t u ) + dm t 


= , (x;y, (Yfy, (zf)', x°, vf, z*, <) - f(t, u)\dt 

-E'[fz(t,u)(k t y + f x (t,u)kt + fy(t,u)(m t )' + f y (t,u)m t + /*(*,«)(«*)' + f z (t,u)n t 
+f v (t, u)(v t - ut)]dt - Z e t dW t . 
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Denote by X^ e = X? + X6(X? + k t ), Y t x ' e = Y t u + X0(Y t e + m t ), Z^ e = + X6(Z? + n t ) and 
u x ' e (t) — u t + X9(v t — u t ). For convenience, we introduce the notation 

/(A) = /(^,(x^ e )^(y ^ ^)^(z^ e )^x^^y ^ A '^z^^^(^)). 

Then, we have 

- e {f{t, (x!y, (Yf y, (zfy, xl Yf, zf, u 6 t ) - f(t, «)) <tt 

= f umxfy + (k t y)d\ + f u\)(x° + k t )d\ + f fy(\)((Y t y + ( mt y)dx 

Jo Jo Jo 

+ [ f y (X)(Y t e + m t )dX+ [ f,(\)((Z e t y + (n t y)d\+[ f z (X)(Z e t +n t )dX 
Jo Jo Jo 



+ / f v (X)(v t - u t )dX, 
Jo 



and Yf satisfies 



- dYf = | jT E/(/ s (A)(A?)' + / X (A)A? + / fi (A)(l»)' + / w (A)i? + f~ z (X)(Z°y + /,(A)2*)dA 

+A + £ t + C t + G t ^dt - Z e t dW u (21) 
where we denote 



A t = 


Te'UMA) 

Jo 


-/*(*,«))(**)' +(/ x (A) - 


f x (t,u))k t } dX, 


B t = 


^0 


- £(*,„))(,„,)' +(/ tf (A) - 


- f y (t,u))m t } dX 


c t = 


^0 


-£(*,„))(,»,)'+ (/.(A)- 


f z (t,u))n t } dX, 


G t = 


/ E'{(/„(A) 
Jo 


- fv(t,u))(v t - u t )} dX. 





A t tends to in L 2 (VL x [0,T]) as 9 — > 0. Indeed, since the Lipschitz continuity of 

f x (x,y,z,x,y,z), 

f x (x,y,z,x,y,z) with respect to (x,y, z,x,y, z), there exists a positive constant C, which 
may differ from line to line if not specified, such that: 

|/x(A)-/x(f,«)| <cxep t , \f x (x)-f x (t,u)\<cxep t , 

with 

I3 t = \(Xf + k t )'\ + \(Y t e + m t )'\ + \(Zf + nt)'| + \X° + fc t | + \Y t 6 + m t | + |Zf + n t \ + h - u t \. 
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Then 



IAI 2 = 1/ E'{(f,(\)-Mt,u))(k t y+(f x (\)-f x (t,u))k t }d\f 

Jo 

< Ce 2 {E'[\(k t )f]E'[^ t ] + \k t \ 2 E'[^} 
= C9 2 (E[\k t \ 2 } + \k t \ 2 )E'[(3 2 }. 



So 



\A t \ 2 dt < C0 2 eQT E[\k t \ 2 ]E'[/3 2 ]dt 

< C9 2 E ( [ \k t \ 4 dt] 2 E ( [ E'[$]dt ) , 



which converges to as 9 — > since the expected values are finite. Similar estimations 
for B t ,C t and G t in (}2~Tj) show that these terms also converge to in L 2 (Q x [0,T]). For 
simplicity, we let I t = A t + B t + C t + G t - Using Ito's formula to \Y t e \ 2 and noting that 
assumption (A2), we have 

E\Y t e \ 2 + E^ \Z e s \ 2 ds < CE \Y s e \\E(X e s )+X e s +E(Y s e )+Y s e + E(Z 9 S ) + Z e s + I s \ds 

< CE \Y s \ 2 ds + ^E^ \Z s \ 2 ds + J e , 



with 



J e = ^J t \X e s \ 2 ds + Ej \I s \ 2 ds, 



where C > is a constant and Jg — > as 9 — > . Applying Gronwall's lemma gives the last 
two results of ([20]). □ 

Lemma 13. Under the assumptions (Al) and (A2), for any v(-) G U, the following varia- 
tional inequality holds: 

E j E'(h x (t)(k t )' + h x (t)k t + hy(t)(m t )' + h y (t)m t + h~ z (t){n t )' + h z (t)n t + h v (t)(v(t) - u^dt 
+E( ly (Y u (0))m } +E(g x (X£)kr) > 0. (22) 
where we denote h(t, (X?)' t u )' t u )' t u , Y t a , Z?, u t ) by h(t). 
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Proof. Since u(-) is an optimal control of the problem, then 



e- 1 



j(u(-) + 9(v(.)-u(.)))-J(u(-)) 



> 0. 



(23) 



From the estimate of (120]) . when 9 — > 0, it follows that 
iE[ 7 (F e (0)) - 7(^(0))] E[ 7y (y«(0))rno], 

±{e jT E' p»(t, (Xf )', (If)', (Z°)', X°, Y t e , Zf, u(t) + 9(v(t) - u(t))) - h(t)} , 
E^ E'(h £ (t)(k t )' + hy{t){m t y + h~ z {t){n t )' + h x {t)k t + h y {t)m t + /^(t)^ + h v {t){v{t) - 



Combining the limits above with 
easily. 



and the definition of the cost functional, we derive 

□ 



4.2 Adjoint equation and Maximum principle 

For deriving the maximum principle, we introduce the following adjoint equation correspond- 
ing to mean- field FBSDEs (jlj), which is a mean-field FBSDEs and whose solution is denoted 

by (p(-),g(0,W). 



(24) 



-dpt = E'[h(t, (x?y t u , Ut )( Pt y + b x (t, (x?)?, u t ) Pt 

+a x (t, (X?)?,ut){q t y + a x (t, {X^)' t \u t )q^dt 

+E' (h x (t) + h x (t) - f x (t, u)(Q t )' - f x (t, u)Q t yt - q t dW t , 

dQ t = E'(fy(t,u)(Q t )' + f y (t,u)Q t - h y (t) - hy{t)}dt 

+E'(/ 2 (t,u)(Q t )' + f z (t,u)Q t - h(t) - h z (t))dW t , 
{ PT = g x (X%) - $ X {X%)Q T , Q = -l y (Y u (0)). 

This equation reduces to the standard one, when the coefficients do not depend explicitly 
on u>' of the underlying diffusion. Under the assumptions (Al) and (A2), this is a linear mean- 
field FBSDEs with bounded coefficients. Moreover, due to Lemma 2.2 and Theorem 4.1 in 
[13J, it admits a unique F-adapted solution (Q,p,q) such that 



E 



sup \Q(t)\ 2 + sup \p(t)\ 2 + [ \q{t)\ 2 dt 

0<t<T 0<t<T Jo 



< +oo. 



Next, we define the Hamiltonian function as follows: 

H(t,x,y,z,x,y,z,p,q,Q,v) 



b(t, x, x, v)p + a(t, x, x, v)q - f(t, x, y, z, x, y, z, v)Q 
+h(t,x,y,z,x,y,z,v). 
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The following theorem constitutes the main result of this section. 

Theorem 14. (SMP in Integral Form). Suppose (Al)-(A2) hold. Let «(•) be an optimal 
control of the problem, and (X U (-),Y U (-), Z u (-)) denote the corresponding trajectory. Then, 
for all v G U, there holds 



E / E' [H v (t, X' t , F/, Z[, X t , Y t , Z t ,pt, q t , Qt, u{t)) (v(t) - u{t))]dt > 0, 
Jo 

a.e.,a.s., where (p(-),q(-),Q(-)) is the the solution of adjoint equation (2~4\ ). 
Proof. Applying Ito's formula to k t pt + m t Qt yields 



(25) 



E 



k T p T + m T Q T - k p - m Qo 
= E 



g x (X%)k T + m ory (Y u (0)) 

{p t b v (t, (X?)', X^ u t ) + q t a v (t, (X?)', X^ u t ) - QJ v (t, u)) (v(t) - u(t)) 



dt 



= E [ E' 
Jo 

-E J E'(k t hx(t) + k t h x (t) + m t hy(t) + m t h y (t) + n t h~ z (t) + n t h z (t)^jdt. 
Together with Lemma 13, we derive 

e / E'\( Pt b v (t, (x?y,x? iUt ) + ot<Tv(t, (x?y,x?,ut) - + Kit)) - «(*)) 

= E f E'[H v (t,(X^\(Yn\(Z^\X^Y t u ,Z^p t ,q t ,Q t ,uA(v(t)-u^ 
Jo 

Thus, we come to the conclusion of this theorem. □ 
Remark 15. From [251) . we can get that 

E' [H v (t, X' t , y/, Z' t , X t , Y t , Z t ,p t , q t , Q t , u(t)) (v - u(t))} > 0, (26) 
dtdP-a.e., for any v £ tl. 

4.3 Sufficient conditions for maximum principle 

This section is devoted to establish the sufficient maximum principle (also called verification 
theorem) of the mean-field stochastic control problem. 
We need the following additional assumptions. 



dt 



(A3) The function $ is convex in x. g is convex in x and 7 is convex in y. 
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Theorem 16. (Sufficient Conditions for the Optimality of the control) Assume that the 
conditions (Al)-(A3) are satisfied. Let u(-) G U with state trajectory (X^,Y t u ,Z^) and 
(p(-),g(-),Q(-)) be the solution of Mean-field FBSDE (24\ ). Suppose 

W[H(t, (X?)f, (Y t u )f, (Z?)/, X?, Y t u , Z? iPt , q t , Q t , u(t))} 
= mmE?)'?y?)?,Y t u ,Z?,p u q t ,Q t ,v)] 

hold for all t G [0, T] . Moreover, suppose function H(t, x, y, z, x, y, z, p, q, Q, v) is convex with 
respect to (t,x,y, z,x,y, z,v) . Then u is an optimal control of problem (ffl)-(Q7|). 

Proof. For any v(-) G U, we consider 

J(v(-)) - J(u(-)) = 1 + 11 + 111 



with 



E [ E' 
Jo 



h(t, (x t 7, 07)', (z?y, xi y", zi v (t)) 



-hit, (Y t u y, (z?y, x?, Y t \ z u tl u (t)) 



dt, 



II 



III 



E 



E 



g(XZ)-g(X$) 



7(^(0)) -j(Y u (0)) 



(27) 



Since g is convex, it holds that 

II = E(g(X v T ) - g(X$) > E[g x (X%)(X v T - X?)]. 
Due to 7 is convex on ?/, we have 

in > 7 ,(^(o))Cr(o) - Y u (0)) = -Q (Y v (0) - y u (o)). 

From Q28J) and (J29D , we get 

II + III > E[g x (X^)(X v T - X") + ly (Y u (0))(Y v (0) - F w (0))]. 
Since $ is convex, 

Y% - Y% = $(X£) - $(X£) > & X (X£)(X%. - X$). 

By applying Ito's formula to Q t (Y t v - Y t u ) + p t (X? - X?), we get 

E[Q T (Y" - F«) - QoCno) - F"(0)) +p T (X^ - X") - p(0)(X v (0) - X«(0))] 
= E[Q T (F^ - y r u ) + 7,(F U (0))(^(0) - Y u (0)) + - $ X .(X£)Q T )(X£ - X£ 



(28) 
(29) 
(30) 
(31) 



2 J E'{ Pt (b(t, {xyy,xiv{t)) - b(t, (x?y,x?, u (t))) - Q t (f(t,v) - /(*,«)) 



+q t (a(t, (X/)', X* ; „(t)) - ff(t, (X«)', X", «(*))) df + IV 



(32) 
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with 



IV 



E 



j\ j (fy(t,u)(Q t )' + f y (t,u)Q t - hy(t) - h y (t))(YX - Y t u ) 

+(Mt,u)(Q t y + f z (t,u)Q t - h{t) - h z (t))(z» - ZX) 
- fat, (x?y, x?, u t ){ Pt y + b x (t, (x?y, x?, Ut ) Pt 

+a x (t, (XX)', X?, u t )(q t )' + a x (t, (X?)', X« u t )q t ) {X? - X") 



~(hx(t) + h x (t) - f x (t,u)(Q t )' - f x (t,u)Q^(X v t - X^dt 



-E 



£ E'{H & (t, (x?y, (Y t u y, (z?y, xx, Y t \ zx, Vt , qt , q u u (t))(xx - x } 



U\l 

t ) 



+H x (t, (xxy, (Y t u y, (z?y, xx, Y t u , zx, Pt , ft , Qu «(*))(*? - x, 

+H y (t, (XX)', (Y t u )>, (ZX)', X?, Y t \ Z?,p t , q t , Q t , u(t))(Y t * - Y t 



n 



+H y (t, (xxy, (Y t u y, (zx)', xx, Y t u , zx, Pt , gt , q u u (t))(Y t v - y?) 

+H g (t, (XX)', (Y t u )', (Z?y,X?,Y t u ,Z?,p t ,q t ,Q t ,u(t))(Zl> - z?y 

+H z (t, (xxy, (Y t u y, (zxy ,xx,Y t u ,zx, Pt , qt ,Q t ,u(t))(zt - zx))dt. 



Together with (l27j). (I30j) . (I3T]) and (132]), we get 
J(v(-)) - J(u(-)) = 1 + 11 + 111 
> I + E J E'{p t (b(t, (XX)', XI v(t)) - b(t, (XX)', XX, u(t))) - Q t (f(t, v) - f(t, u)) 

+q t (a(t, (XX)', XX, v(t)) - a(t, (XX)', XX, u(t)))}dt + IV 
H(t, (XX)', (YX)', (ZX)', XX, Y t \ ZX, Pt , q t , Qt, v(t)) 



E I E' 
Jo 



H(t, (XX)', (1?)', (ZX)', XX, Y t u , ZX, Pt , q t , Q t , u(t)) 



dt + IV. 



(33) 



Noticing H is convex with respect to (x, y, z, x, y, z), the use of the Clark generalized gradient 
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of H, evaluated at ((*«)', (Y?)', (Z?)', X« Z t u ), yields 



#(t, (a?)', (>?)', (2?)', 17, 2?, ft) ft , g t , v(t)) 



-H(t, (X?)', (Y t u ) 


, (z?y 


-yu 


Y t u ,Z?,p t ,q t ,Q t , 


«(*)) 








(z?y, 


X^Y t u ,Z^p tl q tl Q t ,n 


(*))(*?" 


x»y 




)', (z?: 


',x? 


,Y t u ,Z?,p t ,q t ,Qt 


u(ty 




-xt) 




y, (z?: 


',x? 


Y t 1 Z",p t , q t , Q t 


u(t)] 


(Y t v 


_ yuy 


+H y (t, (X?)', Q? 


)', (z?: 


',x? 


Y?,Z u t ,p t ,q t ,Qt 


u(t)] 


(V 




+H s (t, (x?y, ixr, 






Y t , Zt,p t , q t , Q t 


u(t)) 


{z v t 


- z u t y 


+H z (t, (X?)', (Y t u ] 


(z?: 


',x? 


Y t , Zf,p t , q t , Q t 


u(t)) 


{z v t 


-z?) 


+H v (t,(X?y,(Y t u ] 






Y t u ,Z?,p t ,q t ,Q t 


u(ty 


iy t - 


- u t ). 



Combined the above inequality with (1331 . we have 



J{v{-))-J{u{-))> 

E I E' 
Jo 



H v (t, (X?)', (Y t u y, (Z?y, X?, Y t u , Z u t ,p t , q t , Q u u{t)) (v{t) - u(t)) 
Hence, we draw the desired conclusion. 



dt > 0. 



□ 



5 Stochastic maximum principle for fully coupled forward- 
backward stochastic control systems of mean-field type 

In this section, we extend control problems to the fully coupled mean-field FBSDEs. For 
any v(-) G U, the state equation consists of the following forward-backward control system 
of mean-field type: 



dX t = E'[b(t, X' v Y{, Z' t , X t , Y t , Z t , v t )]dt + E'[a(t, X[, Y t ', Z[, X t , Y t , Z u v t )]dW t , 
X(0) = x , 

-dY t = E'[f(t, X' t , Yl, Z[, X t , Y t , Z t , v t )}dt - ZtdWt, 

y T = $(x T ), 



(34) 



where 



b : [0,T] x 
a : [0,T] x 
/ : [0,T] X R x 
$ : R -> R. 



XEX. 

xRx 



x 



XKX. 

! xRx 



X K X 



X 
X 



i d xU ^ R d , 
l d x U ->■ K, 



The expected cost function is given by: 

J(v(-)) = E ( f T E'[h(t, X[, Y^ Z' t , X t , Y t , Z t , v{t))\dt + g{Xx) + 7(Y(0)) 



(35) 
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where 



R R, 
R ->■ R, 
[0,T]x 



x 



X 



X 



x U -> 



The optimal control problem is to minimize the functional </(■) over U. A control that 
solves this problem is called optimal. 
We assume: 

(A4) 

/ (i) b, a, f, $, h, g and 7 are continuously differentiable; 

(ii) The derivatives of b, a, f and $ are bounded; 

(iii) The derivatives of h are bounded by C(l + \x\ + \y\ + \z\ + \x\ + |y| + \z\); 

(iv) The derivatives of g and7 with respect to x and y are bounded by C(l + 
and C(l + |y|) respectively; 

(v) For any given admissible control v(-), the equation (134)) satisfies (H4) and (H5). 

According to Theorem 3.1, for any given admissible control u(-) G tt, there exists a unique 
adapted solution (X%, Y t v , Z%) satisfying the fully coupled mean-field FBSDEs fl34|) . 

Let u(-) be an optimal control and (X u (-), Y u (-), Z u (-)) be the corresponding state trajec- 
tory of stochastic control system. In this case, the corresponding adjoint equation becomes 

' -dp(t) = E'(b £ (t)(pt)' + b x (t)p t + a s (t)(q t )' + a x (t)q t yt 

+E'(h x (t) + h x (t) - Ut)(Q t )' - f x (t)Qt)dt - q t dW t , 

dQ t = E'(fy(t)(Q t y + fy(t)Q t - by(t)(p t )' - by(t) Pt - ~ ^{^t ~ ky(t) ~ hy^dt 

+E'(/ 2 (t)(Q t )' + fz(t)Q t - h(t)(p t y - b z (t)p t - a~ z {t){q t )' - a z (t)q t - h~ z {t) - h z (tj)dW t , 
p T = g x {X%) - ® X (X%)Q T , Q = - ly (Y u (0)), 

in which we use the notation Tp{t) = tp(t, (X^)^")^)^", Y", Zf) for if) = b,a,f,h. When the 
coefficients b, a and / do not depend explicitly on u', the adjoint equation (136|) reduces to 
the standard adjoint equation (see Shi and Wu |21j ) corresponding to fully coupled FBSDE. 

On the other hand, from the assumption (A4) and the fact that (l34p satisfies (H4) and 
(H5), we can easily verify that this adjoint equation (1361) satisfies (H4) and (H6). Then, from 
Theorem 3.2, we know that ( |36|) has a unique F-adapted solution (Q,p,q) such that 



E 



sup \Q(t)\ 2 + sup |p(t)| 5 

0<t<T 0<t<T 



r 



\q(t)\ 2 dt 



< 



-00. 



Define the Hamiltonian function as 

H(t,x,y,z,x,y,z,p,q,Q,v) = 



b(t, x, y, z, x, y, z, v)p + a(t, x, y, z, x, y, z, v)q 

-fit, x, y, z, x } y, z, v)Q + h{t, x, y, z, x, y, z, v). (37) 
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The following two theorems, whose proof are similar to Theorem 14 and Theorem 16, 
respectively and thus are omited, are the main contribution of this section. 

Theorem 17. (SMP in Integral Form). Under assumptions (A4), if u(-) is an optimal 
control with state trajectory (X u (■) ,Y U (■) , Z u (■)) , then there exists a pair (p(-),q(-),Q(-)) of 
adapted processes which satisfies j[3o)) , such that 

E [ E'[H v (t,X' t X,Z' t ,X t ,Y t ,Z t ,p t ,q t ,Q t ,u{t)){v(t)-u{t))}dt>0, (38) 
J o 

P-a.s., for allte [0,T]. 

Theorem 18. (Sufficient Conditions for the Optimality of the Control) Assume the condition 
(A4) is satisfied and let «(•) G U with state trajectory {Xf ,Y t u , Zf) be given such that there 
exist solutions (p(-),q(-),Q(-)) to the adjoint equation j[3h}) . Moreover, suppose the functions 
g, 7, $ are convex and H(t, x, y, z, x, y, z, p, q, Q, v) is convex in (x, y, z, x, y, z, v) . Then, if 

E f [H(t, {x?y, (Y t u y, (z?y, x^ y«, z?, Pt , ftj q u u (t))} 

= mm E'nrynr , y t u , z?, Pt , gt , q u V )\, 



for all t G [0,T], P-a.s., u is an optimal control of problem 



6 Maximum principle for mean-field stochastic games of 
FBSDEs 

In this section, we consider a class of non-zero sum differential games where state variables 
are described by the system of FBSDEs of mean-field type. Our objective is to derive 
necessary conditions for optimality in the form of a stochastic maximum principle and the 
corresponding verification theorem. 

We always use the subscript 1 (respectively, subscript 2) to characterize the variables 
corresponding to Player 1 (respectively, Player 2). 

Let action space C/j be a non-empty, closed and convex subset of IR fc (z = 1, 2, k G N + ). 
The admissible control set is defined as 

U i = {v i eL%{0,T;R k )\v i eU i ,te[0,T}} (2 = 1,2). 

For any Vi(-) £lAi (i = 1, 2), we consider the following mean-field FBSDE: 



dX t = E'[b(t, X' v X t , v 1 (t),v 2 (t))}dt + E'[a(t, X' t , X t , v 1 {t),v 2 {t))]dW u 
X(0) = x , 

-dY t = E'[f(t, X' t , Y t ', Z' t , X t , Y t , Z t , v 1 {t),v 2 (t))]dt - Z t dW t , 

y T = $(x T ), 



(39) 
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where 



b : [0, T] x R x R x U x x U 2 ->■ R, 
u:[0,T]xlxlx[/ 1 xf/ 2 -)' R d , 
/ : [0,T] x R x I x l d x 1 x K x l d x ^ x [/ 2 4 1, 
$ : R R. 

Ensuring to achieve the goal $(xt). Player i (i = 1,2), who has his own benefits, aims at 
minimizing the following expected cost functionals: 

4M0,*2(0) = E [ T E'[h i (t,X^X,Z' v X t ,Y t ,Z t ,v 1 (t),v 2 (t))]dt 

Jo 

+E(^(X T ) + 7i (y(0))), (40) 

where 

& : R ->• R (i = 1,2), 
7i :R-^R(i = l,2), 

^ : [0, T] x R x R x R d x R x R x R d x U x x C/ 2 R {% = 1, 2). 

Suppose each player hopes to minimize her/his cost functional Ji(v\(-), t>2(-)) by selecting 
an appropriate admissible control Vi(-) (i = 1,2). The problem is then to find a pair of 
admissible controls (u\(-), U2{-)) G Wi x W2, called a Nash equilibrium point for the non-zero 
sum game, such that 



■h 
■h 



i(wi(0»«a(0) = V 1 ^ 1 }, J i( u i(0> u 2(0) 

(mi(-), w 2 (-)) = min J 2 (wi(-), w 2 (-)) 
f2(-)eW2 



We call the problem above a forward-backward non-zero sum stochastic differential game 
of mean-field type, where the word "forward-backward" means that the game system is 
described by a FBSDE and the reason for calling "mean-field" is the coefficients of the state 
equation and cost functionals depend on the law of the state process. For simplicity, we 
denote it by Problem (FBNM). 

We assume that the following hypothesis holds. 

(A5) (i) The given functions b(t, x, x, v±, v 2 ), u(t, x, x, v±, v 2 ), f(t, x, y, z, x, y, z, vi, v 2 ), $(x), 
hi(t, x, y, z, x, y, z, vi, v 2 ), gi{x) and j^y) (i = 1, 2) are continuously differentiable with 
respect to all of the components in these functions. 

(ii) All the derivatives in (i) are Lipschitz continuous and bounded. 

For any admissible controls fi(-) and v 2 (-), we suppose that (A5) hold. Then we know 
mean-field FBSDE (J32J admits a unique solution (x Vl > V2 (-), y vl ' V2 (-), z Vl ' V2 (-)) by Lemma 2.2 
and Lemma 2.3, which is called the corresponding trajectory. 
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6.1 A Pontryagin's stochastic maximum principle 

Let (u 1 (-),u 2 {-)) be a Nash equilibrium point of Problem (FBNM) and (X(-),Y(-), Z{-)) be 
the corresponding state trajectory of game system. For any given Vi(-) (zUi (i = 1, 2), since 
Ui is convex, then u\ (•) = + 6{vi{-) - «*(•)) 6 Ui (i = 1, 2), V 9 G [0, 1]. 
We introduce the short-hand notation which will be in force in this section 

6(t) = b(t,X' t ,X t ,u 1 (t),u 2 (t)), 

a(t) = a(t,X' t ,X t , Ul (t),u 2 (t)), 

f(t) = f(t,Xt,Y;,Z' t ,X t ,Y t ,Z t , Ul {t),u 2 (t)), 

hi(t) = h t (t,X^Y^Z' t ,X t ,Y t ,Z t , Ul (t),u 2 (t)). 

Let (k(-) , m(-) , n(-)) be the solution of the following variational equation which is a linear 
mean-field FBSDE: 



h{t){h)' + b x {t)h + b Vl {t){ Vl {t) - u x {t)) + b V2 (t)(v 2 (t) - u 2 {t)) 



dt 



Vx(t)(kt)' + O x {t)k t + CT Wl (t)Mt) - Mt)) + °v 2 {t){v 2 {t) - U 2 (t)) 



dWt. 



dk t = E' 
+E' 

dm f = -E'[/ g (£)(A; t )' + + Jy{t){m t )' + / y (t)m t + f~ z (t)(n t )' + f z (t)n t 

+f Vl (t)(v 1 (t) - ui{t)) + f V2 (t){v 2 (t) - u 2 (t))]dt + n t dW t , 

The adjoint equation corresponding to state trajectory (X Vl,V2 (-), Y Vl,V2 (-), Z Vl ' V2 (■)) , which 
is a mean- field FBSDE and whose solution is denoted by (Pi 1,V2 (-)j Qi 1,V2 (')i QT' V2 ('))> satisfies 



+E'(Mt) + h ix (t) - h(t)(QT' V2 (t))' - Ut)QT' V2 (t))dt - g r> V2 (t)dw t , 

dQT^it) = ^{fy{t)(QT' V2 it)y + fy(t)Q?' V2 (t) - hyit) - h iy (t))dt (^) 

+w(]~ z {t){Q^ v \t))' + f z (t)QT' V2 (t) - ha(t) - h iz (t))dW u 
I PT' V2 (T) = g x (X T ) - $ x (X T )Q^(t), Q = - 7y (y(0)), 

where hi X denotes the partial derivatives of hi with respect to x. 

By (A5) and Lemma 2.2, we can easily verify that the linear FBSDE of mean-field type 



?7> V2 {-l€' V2 {-lQ?> V2 { 

The Hamiltonian function associated with random variables is defined as follows: 



2|) admits a unique solution (p^ 1 ' 2 (-), <lT' V2 (.')■> Qi~' 2 {'))- 



Hi(t, x, y, z, x, y, z,Pi, q h Qi,v 1} v 2 ) 



Pi(t)b(t, x, x, v u v 2 ) + qi(t)a(t, x, x, v 1 ,v 2 ) 
-f(t, x, y, z, x, y, z, v u v 2 )Qi(t) 
+h i (t,x,y,z,x,y,z,v 1 ,v 2 ). (43) 
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Fix«2(-) (respectively, Wi(-)), to minimize the cost functional Ji(vi(-),u 2 (-)) (respectively, 
«^2(wi(-)> v 2 (-))) subject to fl39|) over Ui (respectively, U 2 ) is an optimal control problem of 
mean-field FBSDEs. Following the idea developed in Section 4, it is not difficult to analyze 
the game problem. Thus, we omit the detailed deduction and only state the main result for 
simplicity. 

Theorem 19. (Stochastic Maximum Principle for SDGs) Suppose (A5) hold. Let (ifi(-), u 2 (-)) 
be a Nash equilibrium point for our stochastic game problem (FBNM), (X(-),Y(-), Z(-)) be 
the corresponding trajectory and (Pi 1,V2 ( 4 ), ^'^(O? QT'^i')) ^ e solution of adjoint equa- 
tion ( fjM ). Then we have 

E [ E'[H lVl (t,X' t ,Y^Z' t ,X t ,Y t ,Z t ,p 1 {t),q 1 {t),Q 1 {t),u l {t),u 2 {t)){v 1 {t)-u^ > 0, 

Jo 

E [ E'[H 2v2 (t,XtX,K,Xt,Y t ,Z t , P2 (t),q 2 {t),Q 2 {t),u 1 {t),U2(t)){v 2 {t) - u 2 {t))]dt > 0, 
Jo 

V (fi,f 2 ) eUiX U 2 , a.e.t G [0,T], a.s., 
where the Hamiltonian function Hi is defined by ( [75] ). 



6.2 Sufficient conditions for maximum principle 

We will establish the sufficient maximum principle (also called verification theorem) of Prob- 
lem (FBNM). 

Theorem 20. (Sufficient Conditions for the equilibrium point of Problem (FBNM)) Let 
(A5) hold and suppose that (ui(-), u 2 {-)) G U\ x U 2 with state trajectory (X t , Y t , Z t ) satisfies: 

E'^^x^r/,^,^,^,^,^'" 2 ^),^'" 2 ^),^'^^)^!^),^^))] 

= minE'^^X;,^,^;^,^,^,^^^),^'" 2 ^),^'" 2 ^)^!,^^))], 

E'[iy 2 (t,^,y/,^,x t ,y t ,z t ,p^' W2 (t),^' W2 (t),g™(t), Wl (t), U2 (t))] 

= minE'^^X;,^,^;^,^,^^^^^),^'" 2 ^),^'" 2 ^^!^),^)], 

V2GU2 

for all t G [0, T], where p" 1,U2 (i), g" 1,U2 (i), Q™ 1,U2 (t) is the solution of adjoint equation HM . 
We further assume that the functions $(x), gi(x), 7j(y) and Hamiltonian function Hi (i = 
1,2) are convex in (x, y, z, x, y, z, Vi, v 2 ). Then, (tti(-), u 2 (-)) is an equilibrium point of prob- 
lem (3WS^>- 

Proof. From the proof of Theorem 16, we affirm that 

JiM-W-))- JiM-),u 2 (-))>o 

holds for any ^i(-) G U±, and 

MM-),M-)) = nun MM-),v 2 {-)) 
holds for any v 2 {-) G U 2 . Hence, we draw the desired conclusion. □ 
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Remark 21. Note that if Eq. |32|) does not include the "forward" part and without the 
influence of u' , then the stochastic game problem and the corresponding conclusion reduce to 
the case introduced by Wang and Yu JH^. 



7 Maximum principle for Mean-field stochastic games of 
fully coupled FBSDEs 

In this section, we study the mean-field games of fully coupled FBSDEs. That is, the state 
equation is characterized by following fully coupled FBSDEs: 

dX t = E'[b(t, X' t , y/, Z[, X t , Y t , Z t , v 1 (t),v 2 (t))}dt + E'[a(t, X' t , Yj, Z' t , X t , Y t , Z t , Vl (t),v 2 (t))}dW t , 
X(0) = x , 

-dY t = E'[f(t, X' t , Y>, Z[, X t , Y t , Z t , Vl (t),v 2 (t))]dt - Z t dW u 

y T = $(x T ), 



where 



[0 i r]xRx8xrxRxlxI d xf/ 1 x(/ 2 4 



[0, T] x E x R x R d x R x R x R d x U x x U 2 -> R d 
[0, T I x R x R x R d x R x R x R d x U x x U 2 -> R, 
R — >■ R. 



For any admissible Vi(-) e Hi (i — 1,2), if conditions (H4) and (H5) hold, the fully cou- 
pled mean-field FBSDE flU has a unique F-adapted solution (X Vl ' V2 (-),Y Vl ' V2 (-), Z Vl ' V2 (-)) 
according to Theorem 3.1 . 

For Player i (i — 1,2), the expected cost functionals is defined as follows: 

JiM-),v 2 (-)) = El E'[h i {t,X' t ,Yt,Z' t ,X t ,Y t ,Z t ,v 1 {t),v 2 {t))]dt 

Jo 

+E( c7i (X T ) + 7l (r(0))), (45) 



where 



9i 

li 
hi 



R -»■ R (i = 1,2), 
R R (i = 1,2), 

[0, T] x R x R x M. d x R x R x M. d x ^ x U 2 -> R, (i = 1,2). 



Each player, having the same goal $(Xr), aims at minimizing her/his cost functional 
Ji{vi{-) , v 2 (-)) by selecting an appropriate admissible control Vi(-) £ (i = 1,2). The 
problem is to find a Nash equilibrium point (ui(-) , u 2 (-)) G Wi x W 2 for the non-zero sum 
game, such that 

Ji(iii(0,«2(0) = min Ji(>i(-),m 2 (-)), 

J 2 (u 1 (-),u 2 (-)) = min J 2 (u 1 (-),v 2 (-)). 

v 2 (-)eu 2 
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For simplicity, we denote the problem above by Problem (CFBNM). 

In order to give the maximum principle, we assume that the following hypothesis holds. 
(A6) 

(i) b, a, f, $, hi, gi and r ) i are continuously differentiable; 

(ii) The derivatives of b, a, f and $ are bounded; 

(iii) The derivatives of hi are bounded by C(l + \x\ + \y\ + \z\ + \x\ + \y\ + \z\); 

(iv) The derivatives of gi and7j with respect to x and y are bounded by C(l + \x\) 
and C(l + \y\) respectively; 

(v) For any given pair of control (vi(-) , v 2 (-)) , equation (jUJ) satisfies (H4) and (H5). 

Let (ui(-),u 2 {-)) be a Nash equilibrium point of Problem (CFBNM) and let (X(-), Y(-),Z(-)) 
be the corresponding trajectory of game system. In this fully coupled case, the adjoint equa- 
tion, different from the case in Section 6, has the form: for i = 1,2, 



-d Pi (t) = mlb s (t)(pi(t))' + b x (t) Pi (t) + a x (t)( qi (t))' + a x (t) qi (t)jdt 

+E'(h ix (t) + h ix (t) - Ut)(Qi(t))' ~ Ut)Qi{t))dt - qi(t)dW t , 
dQi(t) = E'(/^)(Q. 4 (t))' + f v (t)Qi(t) - h{t)(pi{t))' ~ k(t)Pi(t) - °y(t)(qi(t)y (4?) 

-<T y (t) qi (t) -hy{t) -hy(t))dt + E'^f s (t)(Q i (t)y + f^Qiit) -b z (t) Pi (t) 
-h(t)(Pi(t))' - a,(t)( qi (t)y - a z (t) qi (t) - h~ z (t) - h tz (t))dW t , 
{ Pi(T) = g ix {X T ) - a> x (X T )Qi(T), QM = -i iy (Y(0)), 

with if) = ip(t, X' t , Y t ', Z' t , X t , Y u Z u m(t) , u 2 (t)) , for ip = b,a, /, hi, h 2 . 

This is a linear fully coupled mean-field FBSDE with bounded coefficients under as- 
sumption (A6). It is easy to know that adjoint equation (HTj) satisfies (H4) and (H6) since 
condition (A6) and equation (154|) satisfying (H4) and (H5). From Theorem 3.2, this equation 
has a unique F-adapted solution (pi(-), %(■); Qi{')) such that 



E 



T 

sup \Qi(t)\' 2 + sup \pi(t)\ 2 + I \ qi (t)\' 2 dt <+oo, (i = 1,2). 

0<t<T 0<t<T 



We define the Hamiltonian function : [0,r]xMxlxM d xlxlxrxlxM rf xlxR fe xM fe ^ 
R as 

Hi(t, x, y, z, x, y, z,p u q h Q h v 1 ,v 2 ) 
= p i (t)b(t,x,y,z,x,y,z,v 1 ,v 2 ) + q i (t)a(t,x,y,z,x,y,z,v 1 ,v 2 ) (48) 
-f(t, x, y, z, x, y, z, v h v 2 )Qi(t) + h{t, x, y, z, x, y, z, v h v 2 ), (i = 1,2). 

The proof of the maximum principle and verification theorem in this case is practically 
similar to Section 5. Thus we present these theorems without proof. 

Theorem 22. (Stochastic Maximum Principle for SDGs of coupled FBSDEs) Let (A6) 
hold. If (ui(-) , u 2 (-)) is a Nash equilibrium point of Problem (CFBNM) and (X (•) ,Y "(•) , Z '(•)) 
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denotes the corresponding trajectory, then for any (vi,v 2 ) G U\ x U 2 , the following maximum 
principle 



E I E' 
Jo 

E E' 
Jo 



H\ Vl (t, X' t , F/, Z^, Xj, It, Z t ,pi(t), qi(t), Qi(t) , Ui(t) , u 2 (t)) (vi(t) — Ui(t)) 

h 2V2 (t, x' t , y/, z;, x t , y t , z t , P2 (t),q 2 (t),Q 2 (t), Ul (t),u 2 (t)) (v 2 (t) - u 2 (t)) 



dt > 0, 
dt > 0, 



hold a.s. a.e., where (pi{-),qi(-),Qi{-)) (i = 1,2) zs £/ie £/ie solution of the adjoint equation 
( Jty and the Hamiltonian function Hi (i = 1,2) is defined by fijty) . 



Theorem 23. (Sufficient Conditions for the Problem (CFBNM)) Assume that the condition 
(A6) is satisfied. Let (ui(-) , u 2 (-)) G U\ x U 2 and (X t ,Y t , Z t ) be the corresponding state 
trajectory. Suppose (pi(-),qi(-),Qi(-)) (i = 1,2) is the solution of linear mean-field FBSDE 
(f^Tp - Moreover, we assume functions $, (i = 1,2) are convex in x, 7^ (i = 1,2) is 
convex in y and function Hi(t, x, y, z, x, y, z, P i, qi, Qi, vi, v 2 ) (i = 1, 2) is convex with respect 
to (x,y,z,x,y,z,Vi,v 2 ). Then, if 

E'[H l (t,X' t X,Z' t ,X t ,Y t ,Z t ,p l (t),q 1 (t),Q 1 (t),u 1 (t),u 2 (t))} 

= mmE'[H 1 (t,X^X,Z' tl X t ,Y tl Z tlPl (t),q 1 (t),Q 1 (t) 1 v 1 ,u 2 (t))] 

Vl&Jl 

E'[H 2 {t ) X[X,K,X u Y u Z t Mt)A2{t),Q 2 {t),u l {t)Mt))] 

= iffinE / [ J ff 2 (t,X t / ) y/,Z t / ,X t) y t ,^,p 2 (i),g 2 (i),Q 2 (t), Ml (t), U2 )] 
V2&U2 

hold for all t G [0,T], (ui(-) , u 2 (-)) is an equilibrium point of Problem (CFBNM). 



8 Applications: Linear-Quadratic Case 

In this section, we give two LQ examples to illustrate our theoretical results. 

Example 24. For notational simplicity, we consider the following one- dimensional stochas- 
tic control problem. Our aim is to search for the admissible control u(-) minimizing 



J(v(-)) = -E 

subject to the following FBSDE: 

dX t = \A(t)E[X t ] + A(t)X t + B(t)v(t) 
-dY t 



v 2 (t)dt + X^ + Y 2 



(49) 



dt + 



C(t)E[X t ] + C(t)X t + D(t)v(t) 



dW t . 



a(t)E[X t ] + a(t)X t + b(t)E[Y t ) + b(t)Y t + 0{t)E[Z t \ + p(t)Z t + E(t)v(t) 
X = a, Y T = X T , t e [0,T], 



dt - Z t dW t , 
(50) 
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where A(-), A(-), B(-), C(-), C(-),D(-), a(-), a(-), &(')> KO.^CO.^CO are bounded 

and deterministic, and v(t), < t <T takes value in R. 

In this process, the Hamiltonian function is in the form of 



-Q 

1 2 

V 



H(t, x,y, z,x,y, z,p, q,Q,v) = p A{t)x + A(t)x + B(t)v + q C(t)x + C(t)x t + D{t)v 

a(t)x + a{t)x + b{t)y + b{t)y + f3(t)z + (5{t)z + E{t)v 

(51) 

where (p(-),q(-),Q(-)) satisfies 



dQ t = \ b(t)E[Q t ] + b(t)Q t )dt + [P(t)K[Q t ] + f3(t)Q t )dW t 

-dp t = (A{t)E\p t ] + A{t)p t + C(t)E[q t ] + C(*)gt - a(t)E[Q t ] - a(t)Q t )<!t - < h d\Y,. 
Qo = —Yq, Pt = Xt — Qt- 

If u(-) is optimal, then it follows from Theorem 5.1 and 137]) that 

u(t) = Q t E(t) - p t B(t) - q t D(t), tE[0,T\. 



(52) 



Moreover, it is easy to check that candidate optimal control [5^) is really the optimal control 
since the coefficients of Eq / f50j) and cost functional ( f^ j satisfy the assumptions of Theorem 
4.2. 



Example 25. Let us consider the following forward-backward stochastic control system: 



dX t 



b(t)E[X t ) + b{t)X t + A(t)E[Y t ] + A{t)Y t + B{t)E[Z t ] + B(t)Z t + D(t)v(t) 



dt 



+ 



P(t)E[X t ] + P(t)X t - B(t)E[Y t ) - B{t)Y t + C(t)E[Z t ] + C(t)Z t + E(t)v(t) 



-dY t = [a(t)E[X t ] + a(t)X t + b(t)E[Y t ] + b(t)Y t + (5{t)E[Z t ] + p(t)Z t + G(t)v{t) 
X = a, Y T = RX T , te[0,T], 



dW t , 

dt - Z t dW t , 
(53) 



where R > is a constant and v G Lp(0,T; U). For simplicity we also suppose that U — R. 
Functions a(-) > 0, o(-) > 0,i(-) < 0, C(-) < 0, A(-) < 0, C(-) < 0, B(-), b(-), B(-), £)(•), E(-), 
G(-), b(-) and /?(•) are bounded and deterministic. For any given v(-), it is easy to show that 
condition (H4) and monotonic condition (H5) hold. Then from Theorem 1, the fully coupled 
Mean-field FBSDEs |33|) has a unique solution (X(-), Y(-), Z(-)). 
The cost functional is 



-E 



L(t)v 2 (t) dt + E[MX% + NY*}, 



(54) 
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where constants M > 0, N > 0. Function L(-) is deterministic and bounded, and L 1 is 
also bounded. By [37]) . the Hamiltonian function is given by 



H{t,x,y,z,x,y,z,p,q,Q,v) 



p b(t)x + b(t)x + A(t)y + A(t)y + B{t)z + B{t)z + D{t)v 
+q\p(t)x + (3(t)x - B(t)y - B(t)y + C(t)y + C(t)z + E(t)v 
-Q a(t)x + a{t)x + b(t)y + b{t)y + f3(t)z + (5{t)z + G{t)v 

+ 1 -L(t)v 2 . 

According to Theorem 17, if u(-) is optimal, then 

u(t) = -L-\t) (p t D{t) + q t E{t) - Q t G(t)) , < t < T, (55) 
where {p{-),q(-),Q(-)) is the solution of the following fully coupled Mean-field FBSDEs 



dQ t = [b(t)E[Q t ] + b{t)Q t - A(t)E[ Pt ] - A(t)p t + B(t)E[q t ] + B(t)q t )dt 

+ (p(t)E[Q t ] + P(t)Q t - B(t)E[Pt] ~ B{t)p t - C(t)E[q t ] - C{t)q^dW t 

-dp t = (b(t)E[ Pt ] + b(t)p t + ~P(t)E[q t ] + 0(t)q t - a(t)E[Q t ] - a{t)Q^jdt - q t dW u 
{ Q = -2NY , p T = 2M T X T - RQ T , t e [0, T]. 

Similarly, it is easy to verify that the monotonic condition (H6) holds, then from Theorem 
10, FBSDEs ( QSP admits a unique solution (Q(-),p(-),q(-))- 

Moreover, since g(x) = M T x 2 , j(y) = Ny 2 , = Rx are convex and H(t, x, y, z, x, y, z,p, q, Q, v) 

is convex in (x,y,z,x,y,z,v), we can know that the admissible control 153]) which satisfying 
the necessary condition of optimality is really an optimal control. 
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